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Abstract 

The superfluid phase transition of bosons in a two-dimensional (2D) system with disorder and an 
external parabolic potential is studied. The theory is applied to experiments on indirect excitons 
in coupled quantum wells. The random field is allowed to be large compared to the dipole-dipole 
repulsion between excitons. The slope of the external parabolic trap is assumed to change slowly 
enough to apply the local density approximation (LDA) for the superfluid density, which allows 
us to calculate the Kosterlitz-Thouless temperature T c (n(r)) at each local point r of the trap. 
The superfluid phase occurs around the center of the trap (r = 0) with the normal phase outside 
this area. As temperature increases, the superfluid area shrinks and disappears at temperature 
T c (n(r = 0)). Disorder acts to deplete the condensate; the minimal total number of excitons for 
which superfluidity exists increases with disorder at fixed temperature. If the disorder is large 
enough, it can destroy the superfluid entirely. The effect of magnetic field is also calculated for 
the case of indirect excitons. In a strong magnetic field H, the superfluid component decreases, 
primarily due to the change of the exciton effective mass. 

PACS numbers: 71.35.Lk, 73.20.Mf, 73.21.Fg, 71.35.-y 
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I. INTRODUCTION 



Superfluidity in a system of spatially indirect excitons (with spatially separated electrons 
and holes) in coupled quantum wells (CQW) has been predicted by Lozovik and Yudson,— 
and several subsequent theoretical studie a 2 i 3 i 4 i 5 i 6 i 7 i 8 i 9 i 10 i n i 12 i 13 i 14 i 15 have suggested that this 
should be manifested as persistent electric currents, quasi- Josephson phenomena and unusual 
properties in strong magnetic fields. In the past ten years, a number of experimental studies 
have focused on observing these behaviors . 16 i 17 i 18 i 19 i 20 i 21 i 22 i 23 One of the appeals of this system 
is that the electron and hole wavefunctions have very little overlap, so that the excitons can 
have very long lifetime (> 20 /is), and therefore they can be treated as metastable particles 
to which quasiequlibrium statistics apply. Also, when the spatial separation between the 
electrons and holes is large enough, the interactions between the excitons are entirely dipole- 
dipole repulsive. 

In real experiments disorder plays a very important ro le . 16 i 17 i 18 i 19 i 20 i 21 i 22 Although the 
inhomogeneous broadening linewidth of typical GaAs-based samples has been improved 
from around 20 meV to less than 1 meV,— the disorder energy is still not always small 
compared to the exciton-exciton repulsion energy. At a typical exciton density of a few 
times 10 10 cm -2 , the interaction energy of the excitons is several meV. On the other hand, 
the typical disorder energy of 1 meV is low compared to the typical exciton binding energy 
of 5 meV, so that the excitons can be viewed as stable particles. Typical thermal energies 
at liquid helium temperatures are ksT = 0.2 — 2 meV. Previous work 2 ^ has shown that 
in the low-temperature limit, in the case of no trapping potential, the density n s of the 
superfluid component in CQW systems and the temperature of the superfluid transition 
(the Kosterlitz-Thouless temperature T<~) decrease with increasing amplitude of the random 
field (the general description of "dirty" boson problem is presented in Ref. 

In the present paper we consider the phase transition in the system of the excitons 
trapped in a confining potential, in which case true Bose condensation is possible. Exper- 
imentally, an in-plane harmonic potential is possible in GaAs quantum well structures by 
using inhomogeneous stress . 17 i 27 The Bose condensation in this case is similar to that for Bose 
atoms in trap . 28 i 29 i 30 i 31 i 32 i 33 i 34 The slope of the external parabolic trap is assumed to change 
slowly enough to apply the local density approximation (LDA) for the superfluid density 
and Kosterlitz-Thouless temperature at each local point of the trap. This assumption allows 
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us to substitute the exciton density profile in the Thomas-Fermi approximation 34 into the 
expressions for the superfluid density and Kosterlitz-Thouless temperature obtained for the 
system without the confinement, found in Ref. We also analyze the phase transitions 
in strong magnetic field, applying the effective Hamiltonian approach developed in Ref. . 

The paper is organized in the following way. In SeclHlthe local density profile, superfluid 
density and Kosterltz-Thouless transition in the trapped indirect exciton system with disor- 
der are obtained. The temperature dependence of the minimal total number of excitons for 
superfluidity at various levels of disorder is calculated. In Sec. IIHI the local density profile 
and the phase transitions of the trapped dirty system in magnetic field are obtained. In 
Sec. |TV]we present our conclusions . 



II. THE DENSITY PROFILE AND THE LOCAL KOSTERLITZ-THOULESS 
PHASE TRANSITION IN LDA 

We consider a slowly changing parabolic potential with the characteristic length of inho- 



mogeneity I much greater than the excitonic de-Broglie wavelength A = h/ ^J2Mp: 

where M = m e + is the mass of indirect exciton (m e and nth are the electron and 
hole masses, respectively), and p is the chemical potential of the system with dipole-dipole 
repulsion, given in the ladder approximation by2S (here and below, h = 1) 

87m 

» = 7 -2 y () 



2Mlog 
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where e is the charge of an electron, D is the distance between quantum wells, s is the 
spin degeneracy, e is the dielectric constant, n is the two-dimensional (2D) exciton density 
without confinement (we consider the dilute exciton gas np 2 1, where p is the exciton 
radius) . 

Following earlier work,^ one can write the random field felt by an indirect exciton in 
coupled quantum wells, induced by fluctuations in the widths of electron and hole quantum 
wells, in the form V(r e ,r h ) = a e [£i(r e ) - £ 2 ( r e)] + ac h [£ 3 (r h ) - ^( r h)}, where separate pa- 
rameters are introduced for the random fields felt by the electrons and holes which make up 
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the exciton, since in the standard coupled quantum well structure the electron and hole are 
spatially separated in two different wells. Here a e ^ = dE^l/dd e} h, where d e> h is the average 
width of the electron or hole quantum well, and E^l is the lowest energy of the electron and 
hole confined state; £1, £ 2 , £3 and £ 4 are the fluctuation amplitudes of the upper and lower 
interfaces of the two quantum wells. We assume that fluctuations on different interfaces 
are statistically independent, whereas fluctuations of a specific interface are characterized 
by Gaussian correlation function (^j(ri)^j(r 2 )) = (7i%5(r 2 — ri), (£i(r)) = 0, where <fc is 
proportional to the squared amplitude of the ith interface fluctuation.— 

As shown in earlier work,— in the low frequency and long wavelength limits, the effect 
of this type of disorder is to give a single parameter Q which is an imaginary self-energy of 
the exciton due to the disorder. The Green's function of the center of mass of the isolated 
exciton at T = in the random field is given in the momentum-frequency domain within 
the second order Born approximation by^ 

G (0) (P^) = , 7^ y (3) 

where \i is the chemical potential of the system, £q(p) — p 2 /2M is the spectrum of the 
center mass of the exciton in the "clean" system, and the parameter Q is determined by the 
disorder as 

«c(0i +92) - a 2 h (g 3 + g 4 ) fm h \ 2 



Q(p,u)=Q = — — —. M; m e = m h ; a h [ — < a e , 

167T 4 \m e J 

Q(p,u>) =Q= g%±^ M, m h » m e ; a h f^V « a e . (4) 

For the slowly changing confinement potential U(r) = / yr 2 /2 (7 is the curvature of the 
confinement potential applied to the center of mass of the exciton taken from Ref . j^l , and r 
is the distance between the center of mass of the exciton and the center of confinement), as- 
suming the inequality Eq. (jlj holds, the Thomas-Fermi approximation can be applied in the 
thermodynamic equilibrium.^ In the Thomas-Fermi approximation, the chemical potential 
of the system at each local point equals the sum of the coordinate-dependent local chemical 
potential and the external field. Since for the "dirty" system the chemical potential enters 



to the interacting Green's function as V — Q 2 (which represents the chemical potential of 
the "dirty" infinite [unconfined] excitonic system^) , we obtain the Thomas- Fermi equation 
for the local density profile n(r) = no — An(r), (where no = n(r = 0) is the exciton density 
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at the point of the minimum of the parabolic potential) 



yJfj?[no] -Q 2 = ^/i 2 [(n - An(r))] - Q 2 + U{r) 



(5) 



where fi[n(r)] represents the contribution to the chemical potential of the "clean" system 
due to dipole-dipole repulsion at each local point, as found by substituting the coordinate- 
dependent local density n(r) into Eq. (J2J). In the very dilute limit (n p 2 — ► 0, p being the 
exciton radius) for the weak confinement potential (An(r) <C n and fi 2 [n ]An(r) / [(p 2 [n ] — 
Q 2 )no] <C 1) we get the approximate exciton density profile 
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A* 2 N - Q' 
/i 2 [n ] 



n r 2 . 



(6) 



The exciton density profiles at different n and Q at T = are represented in Figs. 1 and 2. 

We apply the local density approximation (LDA) to obtain the local superfluid density 
n s (T,r) and the Kosterlitz-Thouless phase transition temperature T c (r) at each point r of 
the trap. The LDA can be used for a potential which is flat on the characteristic length scale 
/ much greater than the average distance R between vortices of the normal phase, where R 
is estimated as^: 

R 17, T r ^/ 2 ' 



In 



cxp 



In 
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P [\ T-T L , 

This estimation works at temperatures above T c , and diverges at T = T c . So we can use 
the LDA at temperatures sufficiently close to T c , but it breaks down at T c . In the LDA 
framework, the trapped system has at each point in space exactly the same properties as an 
infinite system, using the local density n(r) instead of the average density n for the infinite 
system. Substituting the local density profile Eq. © into the expressions for n s and T c for 
the infinite system 2 - 5 , we obtain for the local superfluid density 

3C(3) T 3 n(r)Q 



n s (r) = n(r) 



— s- 



2tt cj(n(r),Q)M 2Mc 2 (n(r), Q) ' 
and for the local Kosterlitz-Thouless temperature 
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Tc(r) 



1 + 



32 /MT»() 



1/3 



Tin'{r) 



+ 1 



32 / MT e °(r) 
\ 27 \ 7cn'(r) 



1/3 



+ 1-1 



??(r) 
2 i/3 ' 



(9) 



where c 2 (r) = J /i 2 [n(r)] — Q 2 /M is the local velocity of sound and s is the spin degeneracy, 
equal to 4 for excitons in GaAs quantum wells. Here T c ° is an auxiliary quantity, equal to 
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the temperature at which the superfluid density vanishes in the mean-field approximation 
(i.e., n s (T c °) = 0), 

, {27in'(r)c 4 (r)M\ 1/3 , s 

T c ■(?) = — l J * } ■ (io) 



In Eqs. © and (Jini), n'{r) is 



3C(3) 



n{r)=n{r) ' S 2Mc^)- (H) 



For the trapped excitonic system we find the dependence of the critical minimal total 
number N c of particles where the global superfluidity exists on the temperature T at different 
Q. The total number of particles corresponding to fixed Uq is given by 

N = 2 n(r)d\ (12) 
Jo 

where, again, the exciton density profile n(r) = uq — An(r) (An(r) is given by Eq. (JBJ)); 
r = yu(2/(7vV 2 ~ Q 2 )) 1 ^ 2 is the root of the equation n(±r ) = 0. 

Expressing n(r) in terms of N using Eq. (j!2j) . and substituting n(r) in Eq. (|5|) . we 
calculate the critical total number N C (T) of excitons corresponding to Kosterlitz-Thouless 
phase transition (see Fig.3). Our calculations show that at fixed temperature T, it takes a 
greater total number of excitons N c to achieve superfluidity as disorder increases. 



III. TRAPPED INDIRECT MAGNETOEXCITONS IN COUPLED QUANTUM 
WELLS IN RANDOM FIELD 

We now consider indirect excitons in a strong magnetic field perpendicular to the quantum 
wells in the presence of the disorder, extending previous work on a homogeneous gas of 
excitons in a magnetic field.— In this case we neglect the transitions between different Landau 
levels of the magnetoexciton, including transitions caused by scattering from the slowly 
varying spatial confinement potential U(r e ,rh) = ^l( r e+ r h) an d the random field potential. 
We also neglect nondiagonal matrix elements of the Coulomb interaction between a paired 
electron and hole. The region of applicability of these two assumptions is defined by the 
inequalities 37 u c ^> E b , to c 3> J (V^ h ^ , where to c = eH/m e _h, and m e _ h = m e mh/(m e + 
rrih) is the exciton reduced mass in the quantum well plane; Eb is the magnetoexciton binding 
energy in an ideal "pure" system as a function of magnetic field H and the distance between 
electron and hole quantum wells D: Ef, ~ e 2 / ' ern\j^ '/2 at D r H and E b ~ e 2 /eD at 
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D th (th = (eif)" 1 / 2 is the magnetic length) . Here (. . .) denotes an average over 
the fluctuations of the random field. The characteristic length of inhomogeneity is assumed 
to be much greater than the magnetic length: / 3> r#. We consider the characteristic length 
of the random field potential L to be much shorter than the average distance between 
excitons r s ~ 1/ ^frm (L< 1/ \/Tvn, where n is the total exciton density) similar to Ref . Q| . 

It can be shown^ that the effective Hamiltonian H e g of the system of indirect "dirty" 
magnetoexcitons at small momenta is identical to the Hamiltonian of indirect "dirty" exci- 
tons without magnetic field but with magnetic mass m# instead of M = me + mh. In strong 
magnetic fields at D > the exciton magnetic mass is uih ~ D 3 e/(e 2 rj I ) 39 . Therefore, for 
the trapped system in strong magnetic field we can apply the expressions for the exciton lo- 
cal density profile, local superfluid density and local temperature of the Kosterlitz-Thouless 
phase transition for the "dirty" trapped (Sec. |H]) system without magnetic field using ran- 
dom field V e ff, and an effective external field of confinement U e g instead of U e (r) + Uh{v), 
respectively. The effective parameter of disorder Q for the system of magnetoexcitons is 
given by2^ 

Q[p,u) = Q = m H . (13) 

and the effective field of confinement has the form 

U eS (R) = \ f exp f- (R ~ r)2 ) [U e (r) + U h (r)} dr. (14) 
7ir H J \ r H J 

Eq. ()14j) is valid if the characteristic length / of inhomogeneity of the trapping potential 
U(r e , Th) is much greater than the magnetoexciton mean size r exc w r#. In a strong magnetic 
field we obtain [7 eff (R) = U(R) = r yR 2 /2. Therefore, magnetic field does not change the 
effective trapping potential. 

Since magnetic field affects the effective Hamiltonian only by replacing the excitonic 
mass M by effective the magnetic mass runr^ the magnetoexciton local density profile can 
be obtained by substitution the effective disorder parameter Q (Eq. (HHJ)) into Eq. © and 
effective magnetic mass m# instead of M and s = 1 into the expression for /i (Eq. (j2J)) 
and Eq. As we can see from this substitution, the magnetoexciton local density at a 
fixed local point r decreases as magnetic field increases proportional to the magnetic mass 
mil ~ H 2 when D 3> r#. The negative correction to the local exciton density is also an 
increasing function of the interwell separation D because ran ~ D 3 . The local superfluid 
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density and the local Kosterlitz-Thouless phase transition temperature can be obtained by 
replacing M by tuh and substitution of the magnetoexciton local density profile in Eqs. (J8J) , 
©, and (|lip. As we can see from Eq. (|13|) . increasing magnetic field H, corresponding 
to higher Q, results in higher minimal total number N c required to observe superfluidity at 
fixed T. 

IV. DISCUSSION 

As we can see from Eq. (J§J), the critical temperature T c for superfluidity at r = is the 
same as that of an infinite system with n = Uq. Since at r > the local exciton density 
decreases, as we go farther from the minimum, the Kosterlitz-Thouless critical temperature 
and superfluid density will decrease. In general, there will be a superfluid centered at r = 
and a normal phase outside this area; the size of the superfluid region will shrink as T 
increases. 

The disorder is found to deplete the superfluid, making the superfluid region smaller than 
it would be in the absence of disorder, according to Eq. (jUJ) (see Figs. 1 and 2). A "clean" 
system is always superfluid at T = 0, but a "dirty" system will have a quantum phase 
transition to the superfluid state at some critical density (see Fig. 3) even at T = 0. 

In a strong magnetic field H, the correction of the local density is found to increase with 
magnetic field as H 2 and with interwell separation D as D 3 . As magnetic field H increases, 
the minimum total number N c required to observe superfluidity increases at fixed T. The 
critical effective disorder parameter increases with magnetic field as H 2 at T = 0. 
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FIG. 1: Exciton density profile at T = for different uq at the disorder parameter Q = O.SSmeV^; 
interwell distance D = 14ram; 7 = 40eV/cm 2 ; dielectric constant e = 13. The parameters of the 
system are taken from Ref.|27|. 
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FIG. 2: Exciton density profile at T = for different disorder parameters Q (in units of meV) at 
no = 3 x 10 11 cm~ 2 ; interwell distance D = 14nm; 7 = 40ey/cm 2 ; dielectric constant e = 13. The 
parameters of the system are taken from Ref. [Z^ . 
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FIG. 3: Critical total number of particles N C (T) corresponding to the Kosterlitz-Thouless phase 
transition as a function of temperature T[K] for different disorder parameters Q[meV] at interwell 
distance D = 14nm; 7 = 40eV/cm 2 ; dielectric constant e = 13. The parameters of the system are 
taken from Ref. 2?1 ]. 
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